Nonanalyticities of entropy functions of finite and infinite systems 
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In contrast to the canonical ensemble where thermodynamic functions are smooth for all finite 
system sizes, the microcanonical entropy can show nonanalytic points also for finite systems, even 
if the Hamiltonian is smooth. The relation between finite and infinite system nonanalyticities is 
illustrated by means of a simple classical spin-like model which is exactly solvable for both, finite 
and infinite system sizes, showing a phase transition in the latter case. The microcanonical entropy 
is found to have exactly one nonanalytic point in the interior of its domain. For all finite system 



sizes, this point is located at the same fixed energy value e c 



jumping discontinuously to a 



different value £j. nfinlte in the thermodynamic limit. Remarkably, ef nltc equals the average potential 
energy of the infinite system at the phase transition point. The result, supplemented with results 
on nonanalyticities of the microcanonical entropy for other models, indicates that care is required 
when trying to infer infinite system properties from finite system nonanalyticities. 
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Phase transitions, like the boiling and evaporating of 
water at a certain temperature and pressure, are com- 
mon phenomena both in everyday life and in almost any 
branch of physics. Loosely speaking, a phase transition 
brings about a sudden change of the macroscopic proper- 
ties of a many-particle system while smoothly varying a 
parameter (the temperature or the pressure in the above 
example). The mathematical description of phase tran- 
sitions is conventionally based on (grand)canonical ther- 
modynamic functions, relating their loss of analyticity to 
the occurrence of a phase transition. Such a nonanalytic 
behavior in a (grand) canonical thermodynamic function 
can occur only in the thermodynamic limit in which the 
number of degrees of freedom N of the system goes to 
infinity Q. 

Different statistical ensembles, like the microcanonical 
or the canonical ones, lead in general to different results 
for statistical averages. Only in the thermodynamic limit 
and under suitable conditions on the interparticle inter- 
actions (namely stability and temperedness 0), equiva- 
lent results are obtained in different ensembles, and one 
speaks of ensemble equivalence (see Q for an introduc- 
tion or 4] for an extensive treatise). For a long time, 
interest in statistical physics concentrated almost exclu- 
sively onto the canonical and the grandcanonical ensem- 
bles, in which explicit calculations typically are easier to 
perform. This focus is mirrored by the fact that even 
the above mentioned and widely accepted definition of a 
phase transition makes reference to these ensembles. 

For large classes of many-particle systems, e. g. in solid 
state physics, the thermodynamic limit is an excellent 
approximation, and one can profit from ensemble equiv- 
alence by performing calculations within the ensemble 
which appears to be the most convenient [lj. However, 



several of the current hot topics in physics deal with sys- 
tems of intermediate size, like polymers, nanosystems, or 
biomolecules. Although these systems are large enough 
for a meaningful statistical treatment, ensemble equiva- 
lence cannot be taken for granted, and it is this fact which 
accounts for the renewed interest in ensemble nonequiv- 
alence and finite system effects. 

In the case of ensemble nonequivalence, instead of 
choosing the statistical ensemble according to conve- 
nience, this choice has to reflect the physical situation 
of interest. When a finite system is, say, energetically 
isolated, it has to be treated within the microcanonical 
ensemble, as was done for example when analyzing the 
experiments on sodium clusters reported in [fj. Then, in 
analogy to the definition of a phase transition as a nonan- 
alyticity of a (grand)canonical thermodynamic function, 
it appears natural to investigate the analyticity proper- 
ties of microcanonical thermodynamic functions. Aston- 
ishingly, little is known about the analyticity properties 
of microcanonical entropy functions of finite systems. 

The very recent observation that the microcanonical 
entropy of a finite system is not necessarily analytic, i. e., 
not necessarily infinitely many times differcntiable on its 
entire domain, came as a surprise even to many experts in 
the field: nonanalytic microcanonical entropy functions 
of finite systems have been reported for classical 0, Q 
as well as quantum systems where the latter case 
relies on a suitable but rather unconventional definition 
of the density of states. Note that such nonanalyticities 
occur for perfectly smooth Hamiltonians, so they are not 
introduced artificially. 

Having observed nonanalyticities in the microcanon- 
ical entropy of a finite system, i.e., in a microcanoni- 
cal thermodynamic function, the similarity to the above 
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mentioned definition of a phase transition becomes ob- 
vious. It may not be adequate to think of such finite 
system nonanalyticities as phase transitions in finite sys- 
tems, and we will back up this statement later in this 
Letter. However, in the same way as nonanalyticities of 
(grand)canonical thermodynamic functions mark "points 
of particular interest" (i.e., phase transitions), the non- 
analytic points of the microcanonical entropy deserve 
special attention and, at least in very small systems, 
should be measurable experimentally. 

Of course not only the microcanonical entropy of a fi- 
nite system can have nonanalyticities, but also its infinite 
system's equivalent. In contrast to the finite system ones, 
nonanalyticities of the infinite system entropy are in gen- 
eral related to phase transitions, and a classification of a 
whole zoo of such transitions can be found in . 

The aim of the present Letter is to discuss the relation 
between nonanalyticities of the microcanonical entropy of 
finite and of infinite systems. This will be done mainly 
by discussing a spherical model with mean field-type in- 
teractions in the presence of a kinetic contribution to the 
energy, but the discussion will be supplemented by re- 
sults from other models as well. This kinetic mean-field 
spherical model is one of the rare examples (the only?) 
which shows a proper phase transition in the microcanon- 
ical as well as in the canonical ensemble and which, in 
both ensembles, allows for an exact solution not only in 
the thermodynamic limit, but also for all finite system 
sizes. These properties make it an ideal laboratory for 
our study of analyticity properties, and to this purpose 
we start by presenting its microcanonical solution. 

Microcanonical solution of the kinetic mean-field 
spherical model. Introduced by T. H. Berlin and M. 
Kac in 1952, the spherical model of a ferromagnet is 
devised such as to mimic some features of the Ising model 
while, at the same time, being exactly solvable in the 
thermodynamic limit for arbitrary spatial dimensions of 
the lattice. We consider a mean field-like simplification of 
the original model where, instead of nearest-neighbor in- 
teractions, all degrees of freedom interact with each other 
at equal strength. The potential energy of this model is 
given by 
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where the N ^ 2 degrees of freedom cr, G M (i = 1, N) 
are subject to the constraint 



N 



= TV. 



(2) 



This constraint restricts the space A at of allowed con- 
figurations a = [ax, . . . , (Tjv) € Aat to an (Af-l)-sphere 



with radius VN. Different from the discrete Ising spin 
variables, the spherical model has a continuous configu- 
ration space. This allows us to endow the system with a 
Lagrangian or Hamiltonian dynamics. Then the kinetic 
energy is given by the quadratic form 



1 N 

Tn{6-) = - X^ 2 ' 



(3) 



where the velocity vector a = . . . ,&n) is an element 
of the tangent bundle of the configuration space Ajy For 
a thus defined model, the microcanonical density of states 
as a function of the energy per degree of freedom e can 
be written as 



f2_/v(e) = ajv / do - / da 

(4) 

with some normalization constant on- The first Dirac 
^-distribution in 10} accounts for the spherical constraint 
©, the second ensures the velocity to be tangent to the 
configuration space, and the third S foliates phase space 
into shells of constant total energy Tjv + Vzv • The form of 
the integral suggests the use of A^-dimensional spherical 
coordinates for both, the a and the a integrations, and 
an appropriate rotation of the coordinate axes renders, 
in each case, all but one of the angle integrations trivial. 
This yields 



fl N (e)(x / drr^" 1 S (r 2 — N) / dfr 



x / dtf sin^ 2 / d^ sin^" 2 



x S (rr cos i?) 6 (r 2 cos 2 d - r 2 + 2Ne) , (5) 

where the trivial angle integrations have been executed. 
Three of the remaining four integrations can be per- 
formed by making use of the properties of the Dirac dis- 
tribution, and the microcanonical density of states can 
be written in the form 



Jo 

x(2e + y) {N - 3)/2 e(2e + y), (6) 

where 9 denotes the Heaviside step function. This in- 
tegral can be expressed in terms of Gamma functions T 
and Gauss hypergeometric functions 2F1 
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< e, 



(7) 



or in terms of Legendre functions. In the interior of its 
support f-i oo), the microcanonical density of states 
r2jv(e) has precisely one nonanalytic point at e 



finite 







for all system sizes N ^ 2, and the same of course holds 
true for the microcanonical entropy defined as 



s N (e) = —lnSl N (e); 



(8) 




FIG. 1: The graph of the microcanonical entropy sjv as a 
function of the energy e, plotted for the system sizes N = 4, 
8, 16, 32, 64, 128, and for the infinite system. 

the first L^T^J derivatives of sat are continuous, where 
[•J denotes the largest integer smaller than or equal to 
its argument. So we see that, although a nonanalyticity 
is present in the microcanonical entropy for any finite 
system size, the strength of the nonanalyticity decreases 
with increasing N. 

Thermodynamic limit. — To compare the above de- 
scribed finite system behavior with that of the infinite 
system, the microcanonical entropy is computed in the 
thermodynamic limit, 



■Soo(e) = Ihn SAr(e) = lim 



N 



lnf2jv(e). 



(9) 



Large N asymptotic expansions of the Gauss hyper- 
geometric functions occurring in (JTJ) can be found 
in for 2F1 (•, ^5^-, •) and in making use 



of the identity Vl-2e 2 Fi (±, - 1; 1 + 2e) = 

I' T> (f^if ) 2 ) • To obtain the leading behav- 
ior only, an evaluation of the integral @) by Laplace's 
method fl3| is sufficient, yielding, apart from a physi- 
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cally irrelevant additive constant, 
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(10) 



see Fig. Q] for a plot. A nonanalyticity in Soo(e) occurs 
i, and we obtain e^ nitc £ e ' c nfinite . Only the 



at e 



infinite 



first derivative of s^, is continuous, so a continuous phase 
transition is found to take place at the critical energy 



see Fig.nfor a plot. The nonanalyticity is such that only £ 



infinite 



Nonanalyticities and the configurational entropy. — In 
a calculation similar to that leading to Eqs. Q and 11UL 
the microcanonical average potential energy (v) can be 
determined as a function of the energy e, and we find a 
critical value of the average potential energy ^■ L ,^ 1 c nnnlte = 
(w)(ej, n tc ) = 0. It may seem a mere chance that this 
value of the critical average potential energy coincides 
with the energy (not average potential energy!) of the 
finite system nonanalytic point, i.e., 



\ infinite 



_finitc 



(11) 



but it is not, and in fact it is a remarkable observation. 
The origin of this coincidence can be understood from 
the configurational entropy s con f of the kinetic mean- 
field spherical model as a function of the potential en- 
ergy v. This quantity, by definition, is identical to the 
entropy of the mean-field spherical model without a ki- 
netic energy term as derived in . The crucial feature of 
Sconf (v) is that it has compact support, and there exists 
a largest value w max = of the potential energy accessi- 
ble for the system. This value w ma x then is responsible 
for both, the nonanalyticity at £^ nltc = in the finite 
system and at e|, nfimtc = | in the infinite system: In the 
finite system, when increasing the energy across the value 
£ = Vmax, suddenly no new configurations a become ac- 
cessible, leading to the nonanalyticity in the finite system 
entropy sjv(e) at e^ nltc = 0. In the infinite system, in 
contrast, the nonanalyticity in Soo(e) is a consequence of 
the average potential energy (v) reaching the value w ma x 
from which on on average no new configurations become 
accessible anymore. As reasoned above, this happens at 
an energy g" 1 '™ 1 ' = I. 

Phase transitions in finite and infinite systems? — We 
remarked in the introduction of this Letter that it may 
not be adequate to think of the nonanalyticities observed 
in the microcanonical entropy sat as finite system phase 
transitions. This statement by now should have become 
clear from the results on the kinetic mean-field spheri- 
cal model: The nonanalytic points of the microcanonical 
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entropy of finite systems become weaker with increasing 
system size N (in the sense that the first discontinuity 
arises in a derivative of order ^). Moreover, their locus, 
i.e., the energy e^ nlto , does not even converge towards 
the critical energy ej, nfinlto of the infinite system. This, in 
our opinion, is not what one would expect from a reason- 
able generalization of the concept of a phase transition 
to finite systems. 

Nonanalytic points of the microcanonical entropy in 
other models. — To get a more complete — and at the same 
time more complicated — picture, we want to supplement 
our findings with results on the number and the quality of 
nonanalytic points of the microcanonical entropy in other 
models. To do so, it is helpful to notice the connection 
between nonanalyticities of the microcanonical entropy 
of finite systems and topology changes within the family 
{Euj^gj^ of constant potential energy submanifolds 



= {a G A N | V N (a) = vN} 



(12) 



In recent studies on the relation between phase transi- 
tions and topology changes in configuration space (see 
[Til ] for an up-to-date list of references), the topology 
changes within the family {X^}^^ have been computed 
for several statistical mechanical models. For different 
models like the mean- field </? 4 model 0], the mean- field 
XY model the mean-field fc-trigonometric model 

[l7j . and the spherical model with nearest-neighbor in- 
teractions 0, a large number of topology changes was 
found, increasing unboundedly with the system size N 
and becoming dense on some interval of potential ener- 
gies v in the thermodynamic limit. Each of these topol- 
ogy changes gives rise to a nonanalyticity in the finite 
system microcanonical entropy. Resorting to Morse the- 
ory, it can be argued on general grounds that, again, the 
strength of the nonanalyticities in sjv decreases with in- 
creasing N such that a first discontinuity shows up in 
a derivative of order ^ 01- These observations make 
the idea of defining a finite system analogue of a phase 
transition via the nonanalyticities of the microcanonical 
entropy — an option which we have already dismissed for 
other reasons — appear even more inadequate: transition 
points, increasing unboundedly in number with the sys- 
tem size, do not seem to match the physical intuition of 
what a phase transition in a finite system should be like. 

Summary. — We have obtained an exact expression 
for the microcanonical entropy of the kinetic mean-field 
spherical model for both finite and infinite system sizes. 
The microcanonical entropy as a function of the en- 
ergy e is found to have exactly one nonanalytic point 
in the interior of its domain. For all finite system sizes 
N, this point is located at the same fixed energy value 



-finite 
* c 

-infinite 



= 0, jumping discontinuously to a different value 

-finite 



Ej™" 1 ™ — \ ^ e^ nltc in the thermodynamic limit of in- 
finite system size. Although at different values of the 
energy, the nonanalytic points in both, the finite system 
and the infinite system cases, are argued to be related to 



the same nonanalyticity in the configurational entropy 
Sconf(w) of the system at a potential energy v = 0. 

Two important conclusions can be drawn from these 
peculiar findings: (i) Finite system and infinite system 
nonanalyticities of the microcanonical entropy are not 
unrelated, but their relation is not a straightforward one. 
Care is required when trying to infer infinite system prop- 
erties from finite system nonanalyticities. (ii) A general- 
ization of the definition of a phase transition via the non- 
analyticities of the microcanonical entropy is too naive an 
approach: the number of such nonanalyticities may in- 
cease unboundedly with the system size, whereas their 
strength may decrease. 
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